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High-temperature superconductivity in the copper-oxide ceramics remains an unsolved 
problem because we do not know what the propagating degrees of freedom are in the 
normal state. As a result, we do not know what are the weakly interacting degrees 
of freedom which pair up to form the superconducting condensate. That the electrons 
are not the propagating degrees of freedom in the cuprates is seen most directly from 
experiments that show spectral weight redistributions over all energy scales. Of course, 
the actual propagating degrees of freedom minimize such spectral rearrangements. This 
review focuses on the range of epxerimental consequences such UV-IR mixings have on 
the normal state of the cuprates, such as the pseudogap, mid-infrared band, temperature 
dependence of the Hall number, the superfluid density, and a recent theoretical advance 
which permits the identification of the propagating degrees of freedom in a doped Mott 
insulator. Within this theory, we show how the wide range of phenomena which typify 
the normal state of the cuprates arises including T— linear resistivity. 
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I. INTRODUCTION 

The secret to solving any many-body problem is to 
correctly identify the propagating degrees of freedom. 
Typically the propagating modes cannot be read off by 
inspecting a Hamiltonian but rather are dynamically 



generated through a collective organization of the ele- 
mental fields. In identifying the principle that leads to 
such organization, it helps to know what to throw out. 
In this context, the Barde en-Cooper-Schrieffer (BCS) 



theorv (|Bardeen et all Il957n of superconductivity in or- 



dinary metals is remarkable because they showed that al- 
though the typical interaction energy scale for electrons 
is on the order of electron volts, only the binding inter- 
action within pairs of electrons, typically of O{10~^eV), 
need be included to obtain a quantitative theory of the 
superconducting state. The underlying principle which 
makes this reduction possible is the resilience of the Fermi 
surface to sh ort-range repulsiv e interact ions. As show n 
by Polchin ski(|Polchinskilll992l). Shankar(IShankaiill994l) . 
and others (jBenfatto and Gallavotti Il990l ). all renormal- 
izations from short-range repulsive interactions are to- 
wards the Fermi surface. As a result, such interactions 
can effectively be integrated out leaving behind dressed 
electrons or quasipart i cles, thereby justifying the key 
Landau tenet! Landau . 19561 ) that the low-energy elec- 
tronic excitation spectrum of a metal is identical to that 
of a non-interacting Fermi gas. Pairing is the only wild 
card that destroys this picture. Since pairing instabilities 
abound in metals for any number of re asons, for e x ample . 



Kohn-Luttinger anomalies ( Kohn and Luttingei . 19651) . 



the F ermi surface is pure mathematical fiction! Laughliru . 
2005[ ) at T = 0. Nonetheless, our understanding of su- 



perconductivity in metals would not be possible without 
it. In this sense, superconductivity within the BCS ac- 
count is subservient to the normal metallic state as su- 
perconductivity emerges as the unique interaction-driven 
instability of the underlying electron Fermi surface. 



Essential then to the success of the BCS theory is a 
clean identification of the natural propagating degrees 
of freedom in a system in which the short-range repul- 
sive interactions may be of arbitrary strength. How- 
ever, there are a number of experimentally relevant sys- 
tems, most notably the copper-oxide high-temperature 
superconductors, in which such an identification of the 
propagating degrees of freedom in the normal state has 
proved elusive, and as a consequence, the nature of 
the superconducting state remains unresolved. This re- 
view focuses on the experimental and recent theoreti- 
cal advances which serve to elucidate the nature of the 
propagating degrees of freedom in the normal state of 
strongly correlated electron systems. By strongly corre- 
lated we mean systems in which no obvious pri i iciple , 
such as that de l ineate d ( Bcnfa tto and GallavottiL Il990t 
Polchinskii Il992t IShankar.. 19941) for Fermi liquids, gov- 
erns the renormalization of the electron-electron inter- 
actions. A ubiquity in such systems in which electron 
absorption is the experimental probe is spectral weight 
transfer over large energy scales, a phenomenon absent 
in Fermi liquids. The presence of UV-IR mixing is the 
tell-tale sign that electrons do not reside in electronic 
states with well-defined energies. Equivalently, the true 
propagating modes are some admixture of multi-particle 
bare electron states. Precisely what is the nature of the 
particles whose energies are sharp is the central question 
in strongly correlated electron physics. More precisely, 
what are the particles for whom the single-particle Green 
function has poles with a non-zero residue? Or equiva- 
lently, what is the natural low-energy theory of a strongly 
correlated electron system? We use natural here to de- 
note a theory in which there are no relevant perturba- 
tions. Knowledge of such constituents would permit a 
straightforward application of the BCS program because 
they render the strongly correlated problem weakly in- 
teracting. The answer to this question ultimately resides 
in the physics of collective phenomena. While collec- 
tive phenomena can arise from states of matter in which 
symmetries are broken, this need not be the case. Noted 
examples include the Kondo effect in which all the elec- 
trons in a metal act collectively to screen a local magnetic 
impurity. Perhaps the example which bares the closet re- 
semblance to the physics here is that of quantum chromo- 
dynamics (QCD). In QCD, the pole in the single-quark 
propagator vanishes at low energy and only bound quark 
states survive. Precisely how such bound states are re- 
lated to the degrees of freedom in the UV is the hard 
problem of QCD. Our key message here is that similar 
physics holds for doped Mott insulators. Namely, com- 
posite excitations emerge as the propagating degrees of 
freedom in doped Mott insulators. In the electron coordi- 
nates, these composite excitations are able to explain the 
UV-IR mixing that is the fingerprint of the strong corre- 
lations in doped Mott insulators as well as the anomalous 
transport in the normal state of the copper-oxide super- 
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FIG. 1 A half-filled band as envisioned by Mott. Each blue 
circle represents a neutral atom with electrons and ground- 
state energy . The energy differences for electron removal 
and addition are explicitly shown. Mott reasoned that no 
doubly occupied sites exist because at zero temperature, U — 
+ ^JV-i _ 2^JV ^ Q ^jjjg course, not true. As 

a consequence the Mott gap must be thought of dynamically 
rather than statically. 



conductors. 



II. MOTT'S PROBLEM: INSULATING STATE OF NIO 



The orig inal Mott dMottl . Il949l ) problem grew out of 
coming to terms with why NiO insulates. Since NiO has 
two half-filled d-levels, it is expected, based on the band 
picture of metals, to conduct at zero temperature. How- 
ever, it insulates. For electrons to conduct, they must 
hop from atom to atom. An impediment to transport 
obtains if the electron repulsions win out. In Mott's con- 
struction, the relevant interaction that dominates for nar- 
row d-bands is the energy cost. 



U 
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E 



N-l 



2E 



N 



(1) 



for placing two electrons on the same Ni atom. Here E'^ 
is the ground state energy for an atom with N valence 
electrons. Excluding the filled levels, N = 2 for Ni. At 
zero temperature, Mott reasoned that there is no Ni atom 
with iV±l electrons if U exceeds a critical value, typically 
on the order of the bandwidth. In such a state, all Ni 
atoms have valence of +2, and no conduction obtains as 
illustrated in Fig. ([Ij. On this account, the resultant 
charge gap is the energy cost for doubly occupying the 
same site with spin up or spin down electrons. As a 
result, the simple Hamiltonian, 



i^Hubb = -tYl 9ij4,aCj,a + 4,t4,i'^hi<^ht^ (2) 

i,(T 



first introduced by Hubbard ([Hubbard, Il963h . in which 
electrons hop among a set of lattice sites, but pay an en- 
ergy cost U whenever they doubly occupy the same site, 
is sufficient to describe the transition to the state envi- 
sioned by Mott. In this model, i,j label lattice sites, gij 



is equal to one iff i,j are nearest neighbours, Cia anni- 
hilates an electron with spin a on lattice site i and t is 
the nearest-neighbour hopping matrix element. In light 
of Eq. ([5]) , the simple Mott picture in which no Ni+++ 
or Ni+ ions exist in the ground state, only works when 
the hopping vanishes. This is the atomic limit. In this 
extreme, the eigenstates of Eq. ([2]) are indexed by the 
number of doubly occupied sites. The propagating modes 
are identified by bringing the interaction term, 

Hu = UY^ rizt^^U " y XI ^1<t'/*<t, (3) 



into quadratic form by defining rjia- 
creates the excitations above the gap. 



6, 



(1 



= cl^rii^a- which 
Its complement, 
r) creates excitations strictly on empty 



sites and hence describes particle motion below the gap. 
Consequently, in the atomic limit, the propagating de- 
grees of freedom can be determined straightforwardly 
from the Hamiltonian. 

Beyond the atomic limit, the propagating degrees of 
freedom responsible for the gap are difficult to pinpoint 
primarily because even the lowest eigenstate of the Hub- 
bard model has doubly occupied character. As a result, 
the charge gap in NiO cannot be thought of in the terms 
envisioned by Mott, namely the gap to the first excited 
state that has doubly occupied character. Stated another 
way, the decomposition of the electron operator as a sum 
of rjia- and ^^o- is not canonical. As a result, rj and ^ 
do not diagonalise the hopping term and hence do not 
propagate independently. Precisely what the propagat- 
ing degrees of freedom are that are the efficient cause of 
the Mott gap is not known. Even in the case of one spa- 
tial dimension where the Hubbard model can be solved 
exactly, it is not tractable to write down explicitly the 
band st ructure of the degre es of freedom that become 



gapped (|Gogolin et all 119981 ) at strong coupling The 



persistence of this problem led Laughlin llLaughlinl 119981 ) 
to assert that the Mott problem and all of its associ- 
ated phenomena, such as the lower and upper Hubbard 
bands are entirely fictitious. He opts instead for anti- 
ferromagnetism as the cause of the gap in a half-filled 
band. Indeed, antiferromagnetism and strong coupling 
Mott physics are closely related as illustrated in Fig. 
Electrons on neighbourin g sites localized b y large on-site 
repulsions can exchange! Anderson , 198?! ) their spins if 
the spins are antiparallel. Second order perturbation the- 
ory around the atomic limit is sufficient to establish that 
the energy scale for this process is J oc 0{t'^/U) as illus- 
trated in Fig. ([2]). As such processes lower the energy, 
long-range antiferromagentism is a natural consequence 
of correlation-induced localization of the electrons pro- 
vided the lattice is, of course, bi-partite. 

However, attributing the charge gap in transition 
metal oxides to symmetry-broken states such as an- 
tiferromagnets leaves an explanatory residue. It i s 
that residue that we term Mottness (|Phillipsl . l2006l) . 
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FIG. 2 Mechanism for the generation of the super-exchange 
interaction when the on-site interaction energy, U exceeds 
the hybridisation energy, t. The exchange energy will scale as 
J = At^/U. 



Consider a prototypical 
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syst em VO2 which 

'2007; 'Biermann etai. 



Good cnough. 
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Koethe et al 
Pouget et a U. '197^ 



1971; Havcrko rt et al\. 
IPaquet and Lerqux-Hugonl 



Qazilbash et all 120081 ; iRice et al. 



2005 : 



198C; 



1994 ; IWent zcovitc h et al\ . 119941 ) a transition to an 
insulating state at roughly 340K. In this system, each 
vanadium ion has a valence of +4 giving rise to a half- 
filled d} configuration. Below 340K, the conductivity 
decreases by four or ders of magnitu d e and a charge 



gap of 0.6ey opens! Qazilbash et al\ . I2OO8I ). While 



no hint of magnetic order is detected in this system, 
the vanadium at oms do pair up t o form tilted dimers 
along the c-axis(lBierma nn et al. \ I2OO5I; 'G oodenoughl . 



1971 



1975 



Paauet and Leroux -Hugon. 1980; P ouget et al. 
Wentzcovitch et all. Il994l). Consequen t ly, 



similar question! Koethe aLl. 2009: Rice aL . 1 199 



Wentzcovitch et all Il994l : IZvlbersztejn and MottI Il975 
has arisen in this system as to whether or not the gap is 
due to symmetry breaking resulting in a doubling of the 
unit cell or the correlation pictur e of Mott. This question 
has p ersisted even though MottI Zvlbersztejn and Mot3 . 
19751 ) pointed out that a gap of 0.6eV is beyond any 



energy scale entailed by dimerization of the vanadium 
ions. The optical response of this system probed by 
ellipsometry is particularly useful here in settling this 
question. The key feature shown in Fig. ([3]) is tha t 
lower ing the temperature to 295Kf Qazilbash et al. 



20081 ). an energy scale considerably less than Q.QeV, 
leads to transfer of spectral weight from states in the 
vicinity of the chemical potential (red region) to those 
at considerably high energies (blue region), roughly 6 
eV away and beyond. Such energy scales over which 
the spectral weight is redistributed vastly exceed those 
relevant to dimerization. They are, however, consistent 
with correlation physics on the [/-scale. This is Mottness 
and it persists, as we will see, even upon a transition 
to the superconducting state. In fact, this state of 
affairs obtains even in half- filled bands, for exam- 
ple, V2O3 , whic h are knowndBaldassarre et aL . 2005 : 
Bao et ali Il997l; [Castellani et all Il978l; lEzhov et al. 



19991 ; iQazilbash et all l2008t iThomas et all Il994l) to 

order antiferromagnetically in the insulating state. The 
second panel in Fig. ([3]) shows the analogous optical 
conductivity across the metal-to-insulator transition in 
V2O3. As in the case of VO2, the transition to the 
insulating state also involves transfer of spectral weight 
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FIG. 3 Real part of the optical conductivity above and below 
the temperature for the onset of the Mott insulating state for 
1) VO2 and V2O3. In both VO2 and V2O3, the transition to 
the insulating state is accompanied by a transfer of spectral 
weight from in the vicinity of the chemical potential to states 
as far as 6eV away. This massive reshuffling of the spectral 
weight upon the transition to the Mott state is a ubiquity of 
the Mott transition. Data reprinted from PRB, 77, 115121 
(2008). 



from the chemical potential to states at least 6eV away. 
This mixing of high and low energy scales is a ubiquity 
in Mott systems. Hence, central to the charge gap are 
propagating degrees of freedom that entail the U scale, 
not simply the smaller energy scale associated with 
whatever ordering phenomenon might obtain. 

UV-IR mixing, as evidenced by the spectral 
weight transfer in the optical conductivity, indi- 
cates that any single-electron description of a Mott 
system is moot. In fact, Mott insulators are 
characterized (iDzyalos hinskii j_ 2003 ;_ Essler and Tsveli"l3 . 
2002t iRoschl [2007: Sta nescu et all l2007l ) by a vanishing 
of the single-particle electron Green function along a con- 
nected surface in momentum space. While the volume of 
this zero surface is not directly tied to the parti cle den- 
sity except in the case of p article-hole symmetry (jRosch . 
2007 : Stanescu et al. . \20ffn ). in direct contrast to the sur- 



face enclosed by the divergence of the Green function 
in a Fermi liquid, the zero surface for a Mott insulator, 
nonetheless, has a profound significance. It indicates that 
electrons are not the propagating degrees of freedom that 
give rise to the gapped spectrum. A correct identifica- 
tion of the propagating modes would result in a pole in 
the associated single-particle Green function. An anal- 
ogy to QCD is in order here. At IR energy scales, the 
single-quark propagator vanishes. However, the meson 
or bound quark propagator has a pole. The Mott prob- 
lem amounts to finding the particle whose pole in the 
single-particle propagator leads to the band structure of 
a Mott insulator depicted in Fig. The only alter- 

native is that some sort of dynamically generated com- 
posite or bound state accounts for the gap in the spec- 
trum. That the natural propagating modes responsible 
for the gapped structure of a Mott insulator are compos- 
ite particles or bound states of the elemental excitations 
can be seen from a simple physical argument. Beyond 
the atomic limit of a Mott insulator, double occupancy 
explicitly occurs in the ground state. Double occupancy 
cannot occur without the simultaneous creation of empty 
sites. If the empty sites move freely, then the Mott insu- 
lator is in actuality a conductor. As is evident from Fig. 

mobile double occupancy will also destroy local an- 
tiferromagnetic correlations. Hence, both the magnetic 
and electrical properties of a Mott insulator demand that 
doubly occupied (doublon) a nd empty sit es (holon) form 
bound states. In fact, Mott ( Mott . 19491 ) anticipated as 
much in so far as the insulating properties are concerned. 
Doublon-holon binding has been observed to lower the 
energ y in variational approaches to the ha l f-fille d Mott 



band! Kaplan et al . 19821 : Yokovama et all 2006). In a 



simil ar vein, Castellani and colleagues (|Castellani et al 



I979I) argued that the Mott state is one in which dou- 



ble occupancy is localized but delocalized i n the metal . 
While there is a tendency in the literature (jGros, Il989l) 
to regard doubly occupied and empty-site bound states 
as simply virtual double occupancy, this is a misnomer. 
The nu mber of doubly o c cupie d sites in the Mott state 
is finite! Castellani et al . I1979I : iKaplan et all [l982) m 
contrast to the Bri nkman-Rice mechanism in which it 
vanishes identicallY (|Brinkman and Ricd . I1970I) . Further, 
as can be seen by diagonalising a small Hubbard clus- 
ter, not all doubly occupied sites in a half-filled band 
occur because they mediate J— scale physics. Conse- 
quently, facing up to the Mott problem away from the 
atomic limit requires an explicit mechanism for the lo- 
calization of double occupancy. Such a mechanism un- 
derlies the propagating degrees of freedom which are re- 
sponsible for the gap and ultimately the onset of anti- 
ferromagnetic order. While it is tempting to invert the 
problem and invoke antiferromagnetism as the mecha- 
nism for the localization of double occupancy, this is 
problematic bec ause it leav es Mottness unexplained as 
noted previously! Andersonl [1997, ). We oflFer here an ex- 
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FIG. 4 a.) Half-filled band in the limit of large f/ > f re- 
sults in a gap in the single particle spectrum. The bands 
shown represent the upper and lower Hubbard bands. In an 
A'^-site system, the total number of single-particle states in 
each band at half-filling is A'^. The electron-removal band, the 
photo-emission (PES) and electron -addition bands, inverse- 
photoemission (IPES) bands correspond to an electron mov- 
ing on empty and singly-occupied sites, respectively. In the 
atomic limit the splitting between the bands is U. b.) Evo- 
lution of the single-particle density of states from half-filling 
to the one-hole limit in a doped Mott insulator in the atomic 
limit of the Hubbard model. Removal of an electron results 
in two empty states at low energy as opposed to one in the 
band-insulator limit. The key difi'erence with the Fermi liquid 
is that the total weight spectral weight carried by the lower 
Hubbard band (analogue of the valence band in a Fermi liq- 
uid) is not a constant but a function of the filling. 



plicit construction of the propagating degrees of freedom 
underlying the Mott gap. 



III. DOPED MOTT INSULATORS: BREAKDOWN OF 
FERMI LIQUID THEORY 

The normal state of the cuprates embodies a panoply 
of phenomena that are inconsistent with Fermi liq- 
uid theory. The most vexing are the st range metal 
characterised by T— line ar resistivity (|Ando et al. 



2004al: iKonstantinovic et all l200l[ ) 



as opposed to the 
and 



quadratic depe ndence pred i cted i n Fermi liquid s, 

the pseudogapi Alloul et al . 19891 : Norman et al. . 1998t 
iTimusk and Stattl . 19991 ) in which the single-particle 
density of states is suppressed, although the super- 
conducting gap vanishes. We highlight these features 
here because the phase diagram of the cuprates. Fig. 
dl]), indicates unambiguously that the correct theory of 
the superconducting state must at higher temperatures 
account for a charge vacuum (that is the elec tronic state) 
that is capable of explaining how the onset (jAndo et al\ . 




hole doping 



FIG. 5 Heuristic phase diagram of the copper-oxide super- 
conductors. In the strange metal, the resistivity is a linear 
function of temperature. In the pseudogap the single-particle 
density of states is suppressed without the onset of global 
phase coherence indicative of superconductivity. The dome- 
shape of the superconducting region with an optimal dop- 
ing level of a;opt ~ 0.17 is quantitatively accurate only for 

La2-a;Sra;Cu04. 



2004al : IKonstantinovic et all 120011) of the pseudogap 
at the temperature scale T* leads to a cessation of 
T-linear resistivity. While numerous theories of the 
pseudogap abound (I Andersopi. 1987 : Chakravarty et al. . 



2001: Franz and Tesanovic, 2001: Kivelson et all 11998: 



Randeria et al. . 1992t Ranninger et al . 1995 ). none offer 
a resolution of the T— linear resistivity problem within 
any realistic model of a doped Mott insulator. Part of 
the problem is that a series of ass ociated phenom ena, 



2000) 



for example, incipient diamagnetis mdXu et al 
dicative of incoherent p airing (iFranz and Tesanovic , 



2001 



2005 



Randeria et al 



1992 



1995 ) , electronic in ho mogen eity dAbbamonte et al. 



Pasupathv et al . 2008t Tranquada et al. 



metry 

2002; Simon and Varma 



Ranninger et al. 



Kivelson et~ai]. Il998t iMachida 



Zaanen and Gunnarssonl. Il989f ) , time-reversal sym- 



1989; 



20041 



2002 



breaking(Fauaue et al . 2006t Kaniinski et al. 



Xia et al. 2008), and 



quantum oscillations (|Doiron-Levraud et all 120071 ) in 
the Hall conductivity, possibly associated with the 
emergence of closed electron (not hole) pockets in the 
first Brillouin zone (FBZ), obscure the efficient cause of 
the pseudogap and its continuity with the strange metal. 
What we propose here is that the degrees of freedom 
responsible for dynamical spectral weight transfer are 
directly responsible for the pseudogap and the transition 
to the strange metal. 
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A. Dynamical Spectral Weight Transfer: More than just 
electrons 

In the electronic state or charge vacuum that accounts 
for the normal state of the cuprates, the key assump- 
tion of Fermi liquid theory that the low-energy spectra 
of the interacting and free systems bare a one-to-one cor- 
respondence must break down. More precisely, the inter- 
acting system must contain electronic states at low en- 
ergy that have no co unterpart in t he non-interacting sys- 



tem. One possibilitv ()Moller et all Il992i) is that spectral 



weight transfer between high and low energies mediates 
new electronic states at low energy that have no counter- 
part in the non-interacting system, thereby leading to a 
breakdown of Fermi liquid theory. We show in this sec- 
tion that this is precisely what obtains in the Hubbard 
model. 

To motivate this pathway for the breakdown of Fermi 
liquid theory, we analyse the Oxyg en Is x-ray absorp- 



tion experiments dChen et alV 119911 ) on La2-a;Srj;Cu04 



(LSCO) shown in Fig. ([SJ. In such experiments, an elec- 
tron is promoted from the core Is to an unoccupied level. 
The experimental observable is the fluorescence yield as 
a function of energy as electrons relax back to the valence 
states. The experiments. Fig. (jS]), show that at a; = 0, 
all the available states lie at 530eV. As a function of dop- 
ing, the intensity in the high-energy peak decreases and 
is transferred to states at 528eV. In fact, experimentally, 
the lower peak grows faster than 2x while the upper peak 
decreases faster than \ — x, x the number of holes. In 
La2-2:Sr2:Cu04, the doping level x can be unambiguously 
determined because each Sr atom produces one hole. The 
separation between these two peaks is the optical gap in 
the parent insulating material. 

The redistribution of the spectral weight seen in the ex- 
periments can be explained within the Hubbard model. 
Since the experiments are probing the flourescence yield 
into the available low-energy states, the relevant theo- 
retial quantity is the number of single-particle addition 
states per site at low energy. 



A 

L = I N{uj)duj, 



(4) 



defined as the integral of the single-particle density of 
states {N{uj)) from the chemical potential, ^, to a cutoff 
energy scale, A, demarcating the division between the IR 
and UV scales. In a Fermi liquid, A can be extended to 
infinity as there is no upper band, whereas in a semicon- 
ductor, A should extend only to the top of the valence 
band to count the states available upon the addition of 
holes. To calibrate this quantity, we compare it with 
the number of ways electrons can be added to the empty 
states created by the dopants. Let this quantity be Uh- 
Consider first the case of a Fermi liquid or non-interacting 
system. As illustrated in Fig. (O, the total weight of the 
valence band is 2, that is, there are 2 states per site. The 




528 530 
PHOTON ENERGY (eV) 



532 



FIG. 6 a) Normalized fluorescence dChen et all Il99ll ) yield 
at the oxygen K edge of La2-iSra;Cu04+i . In the undoped 
sample, the only absorption occurs at 530eV, indicated by B. 
Upon doping the intensity at B is transferred to the feature 
at A, located at 528eV. b) Gaussian fits to the absorption fea- 
tures at A and B with the background subtracted. Reprinted 
from Phys. Rev. Lett. 66, 104 (1991). 




FIG. 7 Evolution of the single-particle density of states in a 
Fermi liquid in the valence band as a function of the electron 
filling, n. The total weight of the valence band is a constant 
2, that is, 2 states per site. Doping simply pushes one state 
above the chemical potential. The integral of the density of 
states below the chemical potential is always the filling, n. 



7 



integrated weight of the valence band up to the chem- 
ical potential determines the filling. Consequently, the 
unoccupied part of the spectrum, which determines L, is 
given by L = 2 — n. The number of ways electrons can 
be added to the empty sites is also = 2 — n (see Fig. 
([7])). Consequently, the number of low-energy states per 
electron per spin is identically unity. The key fact on 
which this result hinges is that the total weight of the 
valence band is a constant independent of the electron 
density. 

By contrast, a doped Mott insulator behaves quite dif- 
ferently. At half-filling the chemical potential lies in the 
gap as depicted in Fig. (jH). The sum rule that 2 states 
exist per site applies only to the combined weight of both 
bands. At any finite doping, the weight in the LHB and 
UHB is determined by the density. As a consequence, 
there is no independent sum rule for the occupied and 
empty parts of each band. That is, there is no indepen- 
dent sum rule for L. Consider first the atomic limit. In 
this limit, the total spectral weight of the lower band. 



1 



^LHB 



N ^ 



(5) 



is given by the anticommutator of the operators that cre- 
ate and annihilate singly occcupied sites. Since each hole 
in a half-filled band decreases the double occupancy by 
one, the weight of the UHB is \ — x. Because the total 
weight of the UHB and LHB must be 2, wc find that 
2 — ri4-l — a; = 2orn = \ — x and m^jj^ = 1 + a; 
in the atomic limit. The weights \ x and \ — x also 
determine the total ways electrons can occupy each of 
the bands. Thus, in the atomic limit, electrons alone ex- 
haust the total degrees of freedom of each band. Further, 
since each hole leaves behind an empty site that can be 
occupied by either a spin up or a spin down electron, 
the ele ctron addition spectru m in the LHB has weight 
L = 2a: (|Meinders et aU 119931 ). Hence, the occupied part 
of the LHB and UHB both have identical weights of 1 — a; 
in the atomic limit as depicted in Fig. (|4]). 

Explaining the experiments fully necessitates going be- 
yond the atomic limit. Away from this limit, the total 
spectral weight of the LHB 



ihlhb = ^ + x+ jjY^gij (flfja) + ■■■ 

ija 



1 + X -f- (ip) 



has t/U corrections (|Harris and Langd . 119671 ) which are 
entirely positive. Here fia are related to the origi- 
nal bare fermion operators via a canonical transforma- 
tion that brings the Hubbard model into block diago- 
nal form in which the energy of eac h block is n il. In 



fact, all orders of p erturb ation theory (jEskes et all 11994 ; 



Harris and Lang^ 1967 ) increase the intensity of the 
LHB beyond its atomic limit oil + x. It is these dynam- 
ical corrections that a denotes. That a is positive can 
be seen from the simple fact that turning on the hopping 



increases the total weight of the LHB from the atomic 
limit of one per site to ultimately two per site in the 
non- interacting limit. This increase beyond the atomic 
limit of the intensity of the LHB is significant because 
the number of ways of assigning electrons to the LHB 
remains fixed at 1 -I- a;. As a result, electrons alone do 
not exhaust the degrees of freedom in the LHB, in direct 
contrast to the atomic limit. That is, a low-energy theory 
of the LHB requires new non-fermionic charge degrees of 
freedom. Nonetheless, there is a conserved charge given 
simply by the electron filling. Since the low-energy the- 
ory must have both fermionic and non-fermionic degrees 
of freedom, there are clearly less fermionic quasiparticles 
than there are bare electrons. Consequently, the Lan- 
dau Fermi liquid one-to-one correspondence between the 
two fails. We propose that the chemical potial for the 
effective number of low-energy fermionic degrees of free- 
dom can be determined by partitioning the spectrum in 
the LHB so that dynamical spectral weight transfer is 
essentially removed. In such a picture, the empty part 
of the spectrum per spin is equal to the weight removed 
from the occupied part of the LHB when a hole is cre- 
ated. Hence, we arrive at the assignments of the spectral 
weights in Fig. ([8}d) in which the doping level is renor- 
malized by the dynamics; that is, x' = x + a. In other 
words, the dynamical degrees of freedom denoted by a 
serve to supplement the effective phase space of a hole- 
doped system and x' — x + a now denotes the effective 
number of hole degrees of freedom per spin at low energy. 
Any experiment that couples to the fermionic low-energy 
degrees of freedom should be interpreted in terms of the 
total number of hole degrees of freedom, x + a not x. In 
terms of the true fermionic quasiparticles, L = 2x' which 
should be compared with rih = 2x. Clearly, L/rih > 1 
(as is seen e xperimentally) and Fermi liquid theory fails. 
This failure (|Chakrabortv "oil . l2009f ) arises ultimately 
because the chemical potential for the fermionic degrees 
of freedom which make the doped Mott system weakly in- 
teracting is less than that of the bare electrons. Precisely 
how Fermi liquid theory r e-emerges beyond a critical dop- 



ing is detailed elsewhere (|Chakrabortv et all 12009) 



Experiments on the temperature dependence of the 
Hall coefficient offer direct confirmation that the dop- 
ing level is renormali zed dynamically. In LSCO , the H all 
number has been fit (jGor'kov and TeiteFbauml . 120061 ) to 
a function of the form, 



"Haii(r,a:) = no[x) -|- ni(a;) exp[-A(a:)/r] (7) 



where n njx) is temperatur e independent. 
Empiricallv (|Gor'kov and Teitel'bau^ . l2006l ) A(x) 
is related to the pseudogap scale. Consequently, there 
seem to be two types of charges in LSCO, ones that 
arise from the doping and others which emerge from 
the transfer of anti-bonding states down to the chemical 
potential producing effectively bound charge states. 
Our key contention here is that a arises from the 
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Doped Mott Insulator 



B. Optical Conductivity: Emergence of the pseudogap 
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FIG. 8 Redistribution of spectral weight in the Hubbard 
model upon doping the insulating state with x holes, a is the 
dynamical correction mediated by the doubly occupied sec- 
tor. T o order t/U, this correct ion worked out by Harris and 
Langeplarris and Langel.ll967t). a) The tradit ional approach 
(jHybertsen et all 11993 : iMeinders et all 119931 ) in which the 
occupied part of the lower band is fixed to the electron filling 
1 — X. b) New assignment of the spectral weight in terms of 
dynamically generated charge carriers. In this picture, the 
weight of the empty part of the LHB per spin is the effective 
doping level, x' = x + a. 



second term in Eq. ([7]). In part, it is the presence 
of these two distinct kinds of charge carriers, inferred 
from the Hall number, th at has motivated a two- 
fluid model of the cupra tes( Barzvkin and Pined . l2006t 
Gor'kov and Teitel'baum . 20061) . One of the key points 
of this Colloquium is that two types of charges are 
already implied by L > 22; as noted in Fig. ((8)3). The 
extra degrees of freedom arise from the hybridization 
with the doubly occupied sector. Interestingly, the 
same hybridization arises even in the half-filled band. 
However, a gap appears. Accounting for this difference is 
a charge 2e boson which forms bound states with a hard 
gap at half-filling but only a pseudogap in the doped 
case as implied by Eq. ([T]) . This new degree of freedom 
emerges only when the high-energy sector is integrated 
out exactly. In fact, this procedure proves quite generally 
that integrating out the high energy scale in a doped 
Mott insulator and in a Fermi liquid is fundamentally 
different. In the latter, no new degrees of freedom are 
generate d, whereas in the fo r mer a new charge 2e boson 
emerges Jchov et all l2008allbl: iLeigh and PhilU^ , l2008t 



Leigh et all 120071: IPhillips et aU . l2009t ). The charge 2e 
excitation mediates new electronic states at low energies 
by binding to a hole and hence explains naturally the 
temperature-dependent component of the Hall number. 
Such bound states generate a pseudogap and their 
unbinding yields T— linear resistivity. 



Optical conductivity experiments also indicate that the 
effective number of charge carriers exceeds the nominal 
count provided by the doping as predicted by the quasi- 
particle picture in Fig. ([HJd). To see this, we compute 
the effective number of carriers, or more precisely the 
normalized carrier density, by integrating the optical con- 
ductivity 



2mVr 



cell 



7re^ 



a{u!)dLu 



(8) 



up to the optical gap il ~ 1.2eV. Here a{uj) is the optical 
conductivity, Vcoii the unit-cell volume per formula unit, 
m the free electron mass, and e the electron charge. In a 
rigid-band semiconductor model in which spectral weight 
transfer is absent, i Vpff = nh- Sliqwn in Fig. @ is the 



optical conduc t ivitv(|Cooper et all Il993t iLee et al. ■ l2005t 



Ortolani et all l2005al) for YBa2Cu30y along the Cu02 
plane for oxygen doping levels of y = 6.1 and = 6.6. At 
y — 6.6, the strong suppression of the optical conduc- 
tivity below ss 1.2eV vanishes and is accompanied by a 
decrease in the spectral intensity in the high-energy sec- 
tor. This behaviour is analogous to the inversion of the 
spectral weights above 340K in the VO2 upon a transi- 
tion to the metallic state. In fact, the energy scales in 
the doping-induced metallic state are identical to that 
in the temperature-induced transition in these two sys- 
tems. To quantify the transfer of spectral weight, we 
plot the effective number ()Cooper et "all. ll990l) of carriers 
iVcff (shaded red region) that fill-in the optical gap in the 
'parent' material as a result of doping. As is evident from 
Fig. (fTO|). the normalized carrier density exceeds the car- 
rier density that would obtain in a doped semiconductor 
model (the dashed line). This re sult is significan t be- 
cause independently it was shown! Lee et al . 2005f) that 
throughout the underdoped regime of the cuprates, the 
effective mass is constant. As a result, the Mott tran- 
sition proceeds by a vanishing of the carrier number 
rather th an the mass divergence o f the Brinkman-Rice 
scenario ljBrinkman and Rice . 1970() . Consequently, the 
strong deviation from 2x is a clear indicator that dy- 
namical spectral weight transfer is operative as well in 
the optical conductivity. This is expected since the cur- 
rent couples to the true low-energy degrees of freedom 
and hence the quasiparticles identified in Fig. ^p) which 
have an effective doping level of x' = x + a > x are rele- 
vant. A final feature which deserves mention is the non- 
zero intercept of A^off at a; = 0. This suggests that even 
at half-filling, a remnant of the excitations that fill-in the 
spectral density below the Mott gap is present. 

The precise nature of the excitations that lead to a 
fiUing-in of the spectral weight immediately above the 
lower Hubbard band is intimately tied to the pseudo- 
gap. As such states emerge from high energy, they 
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FIG. 9 Evolution of the in-plane optical conductivity in 
YBa2Cu30y (YBCO) for two doping levels. To a good fac- 
simile, y = 6.1 represents the parent material and y = 6.6, 
an overdoped sample. The key feature is the transfer of spec- 
tral weight above 1.2eV in the 'parent' material to states in 
the vicinity of the chemical potential upon doping. Reprinted 
from PRB, 47, 8233 (1993). 



arise from the incoherent part of the spectrum and 
hence must necessarily be anti-bonding and hence non- 
current carrying. That is, they describe localized exci- 
tations. This is consi stent with the fac t that the oper- 
ators which generate (|Eskes et al. 1, Il994l) the dynamical 
part of the spectral weight are inherently local. Alter- 
natively, L > 2x means that while there are L ways to 
add a particle, there are only 2x ways to add an elec- 
tron. This mismatch means that some states must be 
orthogonal to the addition of an electron. Consequently, 
the spectral function at certain momenta contains states 
in which the spectral weight vanishes arbitrarily close 
to the chemical potential. For L — Uh of the particle 
addition states, a quasiparticle peak is absent in the 
electron addition spectrum. The result is a pseudogap 
as all k-states are not necessar ily gapped. Although it 



has been known for some time (|Harris and Langd . 119671 ) 
that L > 2a; at strong coupling in a doped Mott insula- 
tor based on the Hubbard model, that this simple fact 
implies a pseudogap has not been deduced previously. 
This conclusion is borne out experimentally. While the 
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FIG. 10 Integrated optical conductivity for a hole-doped 
(YBCO) and an electron-doped (Nd2-xCe2Cu04 (NCCO)) 
cuprate reflecting the normalized carrier density. The dashed 
line indicates what is expected for a doping a semiconductor. 
Within the Hubbard model, the charge density in excess of 
the nominal doping level is generated by dynamical spectral 
weight transfer. Reprinted from PRB, 41, 11605 (1990). 



pseudogap is most easily seen through the c-axis opti- 
cal conductivity, it can also be detected, though more 
subtlety, from in-plane transport measurements. A s the 
lower panel in Fig. ((TT]a) indicates (|Lee et aLl . lioOSi) . the 
optical conductivity displays broad features at interme- 
diate energies with a well-defined peak at zero frequency. 
This peak is well-described by a dependence in a 

Drude model. However, below T*, a distinctive band 



appears (|Lee et all 120051 ) in the optical conductivity in 
the mid-infrared energy range. The central frequency of 
this mid-infrared feature evolves to lower energies and 
gradually disappears in he avily overdoped samples. The 
evolution! Lee et al ]. l2005h of wmir as a function of dop- 



ing (see Fig. (|llb ) tracks well that of T*, thereby cor- 
roborating that the mid-IR feature and the pseudogap 



IR remains a 


point of controversy! Chakrabortv et al. 


2008: Coman 


ac et al. 2008: Haule and Kothar. 2007 


Hwang et all 


2007; Moore et all 20021) in cuorate ohe- 



nomenology, two things are clear. First, any explana- 
tion of it must apply equally to the pseudogap regime. 
Second, the correct explanation must involve dynam- 
ical spectral weight transfer from the upper Hubbard 
band as it is from this band that the mid-IR spectral 
intensity originates. The latter implies that the mid- 
IR should be absent from low-energy reductions of the 
Hubbard model which ignore double occupancy. It is 



for th is reason that a recent analvsis (|Mishchenko et al. 



20081 ) of the mid-IR within the t-J model has found it 
necessary to invoke phonons to account for this reso- 



10 



nance. In fact, Uchida, et al. dUchida et aZj . Il99l[ ) were 
the first to point out that because the mid-IR scales 
as t, such physics is beyond that captured by the t-J 
model, at least in its traditional implementation in which 
the operators are not transformed (see first Appendix). 
However, in the Hubbard model, be it the s ingle- 



(a) 



banddChakrabortv et al. 



, .20081 IComanac ad l2008t 

Ortolani et aLl. l2005bl: iToschi et aLl. l2005l) or t he multi- 
band Hubbard model! Haule and Kotliaii 2007 ). despite 
the differences in the computational scheme, a mid-IR 
resonance appears. Withi n the computational sche me 



used by Chakraborty, et al. (jChakrabortv et aLl . 120081 ). it 



is possible to determine which local correlations on a pla- 
quette contribute to the mid-IR. They found that of the 
256 plaquette states, when just a single state was elimi- 
nated, the mid-IR feature vanished. Key features of this 
state are that 1) it contains a mixture of singly (87%) 
and doubly occupied (13%) sites, 2) its energy is — 1.3t, 
essentially the energy of mid-IR peak, and 3) the spatial 
symmetry of the eigenstate is dx^-y2. The latter is par- 
ticularly revealing because it indicates that the mid-IR is 
highly anisotropic and has the same momentum depen- 
dence as doe s the pseudogap. thereby corroborating the 



experimental (iLee et ali . |2005[ ) trend in Fig. (fTTb ') that 



the pseudogap and mid-IR share the same origin. 



C. Scale Invariance: Beyond one parameter 

Scale invariance is a fundarnental prope rty of critical 
matter. For quantum matter (jHerta . 119761 ) that is crit- 
ical, scale invariance implies that the system looks the 
same on average at any chosen time and length scale. 
Provided that w < T, quantum critical systems exhibit 
classical behaviour in which only the temperature gov- 
erns the dynamics in the vicinity of the critical point. In 
this quantum critical regime, as it is called, the transport 
relaxation rate between the quasiparticles is universal in 
that in can be deduced. 



1 

Ttr 



knT 



(9) 



by using dimensional analysis on Boltzmann's constant, 
T and h, Planck's constant. Coupled with the assump- 
tion of scale invariance, we posit as well a Drude form 
for the optical conductivity. 



fDrudc 



4:TT1+ UJ^Tt 



2_2 ' 



(10) 



with cjpi the plasma frequency. These two assumptions 
lead immediately to T— linear resistivity when Eq. (|9]) 
is substituted into Eq. ([T0| and the zero-frequency limit 
is taken. It is for this reason that T— linear resistiv- 
ity is commonly attributed to quantum criticality. Two 
questions are relevant here, however, to determine if this 
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FIG. 11 a) Optical conductivity for LSCO (left panel) and 
YBCO (right panel) as a function of frequency in the under- 
doped regime. The top two contain measurements of the opti- 
cal conductivity at lOK for nonsuperconducting samples, the 
next four show cr(tj) at T = Tc for superconducting samples, 
whereas in the bottom panel, the temperature is roughly T* 
for the y = 6.75 material. In all spectra, phonons have been 
removed by fitting them with Lorentzian oscillators. Clearly 
shown in the underdoped samples is a Drude-like peak (red) 
followed by a resonance in the mid- infrared (blue). The ap- 
pearance of the mid-IR below T* indicates a strong connec- 
tion with the pseudogap. b) Frequency of mid-infrared, cjmir, 
and the pseudogap onset temperature T* as a function of dop- 
ing. The combined onset of the mid-IR below T* and the fact 
that a;MiR and T* have the same doping dependence indicates 
that they share a common cause. Reprinted from, PRB, 72 , 
054529 (2005). 
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procedure is internally consistent with the underlying hy- 
pothesis of scale invariance: 1) Does Eq. (|10p describe 
the cuprates in any doping range, particularly at optimal 
doping and 2) Is Eq. ([TU]) consistent with scale invariance 
of the conductivity? To answer the first question, we con- 
sult the data shown in Fig. ([T^ in which the optical con- 
ductivity in optimally doped Bi2Sr2Cao.g2Yo.o8Cu208+5 
is plotted assuming the Drude formula is valid. Specif- 
ically, it is assumed that the conductivity scales as 
T-f'fiuj/T) as in the Drude formula with l/rt,- cx T. 
Instead of a universal func tion for the entire frequency 
range, van der Marel, et al. (|Marel et "all. l2003l) find that 
/X = 1 for Lu/T < 1.5 and fi w 0.5 for uj/T > 3 if /(w/T) 
is described by the Drude formula. While such a change 
in the exponent fi is inconsistent with the Drude for- 
mula, there is a deeper point here. Namely, the Drude 
formula is not a general consequence of scale invariance. 
The requirements that quantum criticality places on the 
for m of the conductivity have been explicitly worked 
out (|Phillips and Chamonl . l2005h under three general as- 
sumptions, 1) the charges are critical, 2) there is only 
one critical length scale and 3) charges are conserved. 
Under t hese three general condition s, the conductivity 
acquires! Phillips and Chamonl . 12005 ) the form. 



(11) 



where Q is the charge, S is a universal function, and z is 
the dynamical critical exponent defined as the exponent 
that relates time and space. Technically, the dynamical 
exponent can be defined as 



dlne(p) 
dh\p 



(12) 



where e(p) is the quasiparticle dispersion as a function 
of momentum. If space and time are on equal footing, 
then z — 1. Note absent from Eq. pT|) is any addi- 
tional energy scale such as the plasma frequency. Hence, 
a conductivity of the Drude form is incompatible with 
strict scale invariance within a single-parameter scaling 
scenario. As a consequence, while the cuprates might be 
describable by some kind of quantum critical scenario, it 
is not the naive one in which one-parameter scaling is op- 
erative. The root cause of this is the mixing between the 
high and low energy scales which pervades the normal 
state of the cuprates. 



D. Superconducting State: Color change 

The phenomenon of UV-IR mixing detailed in the pre- 
ceding sections might on the surface be thought to be 
irrelevant once superconductivity obtains. That is, only 
phenomena on the energy scale of the pairing interaction 
should be relevant to the ground state of the supercoduct- 
ing cuprates. However, this is not the case. Riibhausen, 




FIG. 12 Temperature/frequency scaling behaviour of the real 
part of the optical conductivity of Bi2Sr2Cao.92Yo.o8Cu208+a 
(BSCO). The data are plotted as a function of 
(aj/cj!j-^cr(aj,r)). Using a function of the form, T-''f{u}/T) , 
we observe a) collapse of all curves for fj, — 0.5 for huj > ksT 
and b) collapse for /i = 1 for fiu/kBT < 1.5. Such a change in 
the exponent ^ is not consistent with the general assumption 
of scale invariance. Reprinted from Science, 425, 271 (2003). 



et al. (|Riibhausen et all l200l[ ) observed, that in under- 
doped BSCO, changes occur in the optical conductivity 
up to 3eV or lOOA, where A is the superconducting order 
gap. But perhaps the changes in the 3eV range are just 
indicative of some strong-coupling effect that has noth- 
ing to do with the condensation to the superconducting 
state. To answer this question, we focus on the f-sum 
rule, 

A = 7meV(2m) = / a{uj)duj. (13) 
Jo 

In understanding the spectral changes in a high-T^ super- 
conductor, it is helpful to separate A into a low-energy 
component. 



Al — I a[uj)duj 
Jo+ 



and a high-energy part 

Ah = 



a{uj)duj. 



(14) 



(15) 



The cutoff is chosen so that Ah contains strictly the 
spectral weight associated with interband transitions. 
Typically, r2/(27rc) — 10, OOOcm"^ is sufficient to demar- 
cate the minimum of cr(w) which demarcates the bound- 
ary between the intra-band and inter-band transitions. 
The opening of a gap opens in the optical conductivity 
accompanies the transition to the superconducting state. 
The spectral weight removed for fioj < A, A the super- 
conducting gap, is transferred to a (5— function at zero 
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frequency. The weight in t he 5 funct i on is captured by 
the Ferrell-Glover-Tinkham( Tinkhaml 120041 ) sum rule 



D = Al^- 



4* 



(16) 



In BCS superconductors, there is no contribution to D 
from Ah- Typically, D is recovered simply by inte- 
grating up to no mo re than lOA. However , the ellip- 



changes in the dielectric function obtain for lOOA sug- 
gest other wise for the cuprates . Indeed this is so. 



Santander-Svro et all 120031 ) 



For BSCO (|Molegraaf et aZ.1 . 120021 : lOrtolani et aLl . l2005bD 
both optimally and underdoped, diminishes as the 
temperature decreases and a compensating increase is ob- 
served for Ai as depicted in Fig. ([13]). This indicates that 
it is the loss of spectral weight in the high-energy sector 
that drives the superconducting state. These data are 
also consistent with other optical measurements which 
indicate that the full weight of the 5— function in the su- 
perconducting state is reco vered by integrating the opti - 
cal conductivity out to 2eV( 
and numerical calculations 
frequency-dependent pairing interaction in the Hubbard 
model involves a non-retarded part that arises entirely 
from the upper Hubbard band. This color change from 
the visible to the infrared implies that superconductiv- 
ity in the cuprates is fundamentally different from that in 
metals. That is, in the cuprates, superconductivity is not 
simply about low-energy physics on a Fermi surface. The 
correct theory should explain precisely how loss of spec- 
tral weight at high energies (2 eV away from the chemical 
potential) leads to a growth of the superfluid density. 
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FIG. 13 Temperature dependence of the low-frequency spec- 
tral weight Ai+_D(r) and the high-frequency spectral weight 
ylh(r) for optimally doped (top) and underdoped (bottom) 
Bi2Sr2CaCu208-5. The insets show the derivatives of these 
quantities multiplied by . Reprinted from Science, 295, 
5563 (2002). 



A. Traditional Method 



IV. WILSONIAN PROGRAM FOR A DOPED MOTT 
INSULATOR 

The essence of the optical experiments on the normal 
state of the cuprates is that the number of particle addi- 
tion states per electron per spin exceeds unity, in direct 
violation of the key Fermi liquid tenet. Within the Hub- 
bard model, this state of affairs obtains because of the 
dynamical mixing of the UHB and the LHB. That is, it is 
absent if double occupancy of bare electrons is projected 
out. In this limit, L = 2x and L/rih — 1. The key ques- 
tion that arises is: How can such mixing be incorporated 
into a low-energy theory of a doped Mott i nsulator? As 



has been pointed out prev i ously, L/uh > idLeigh et i 
20071 : iMeinders efoZI . Il993l : IPhillips oil . l2009|), or more 



generally that the intensity of the LHB exceeds 1 -I- x, 
implies that the true low-energy theory of a doped Mott 
insulator must contain more than just electrons. As the 
new degrees of freedom arise from the mixing with doubly 
occupied sites, it is reasonable that some sort of collective 
charge 2e excitation should emerge. We show explicitly 
here how this physics arises. 



Of course, the traditional approach to constructing 
a low-energy theory of a doped Mott insulator does 
not involve isolating the new degree of freedom respon- 
sible for dynamical spectral weight transfer. Rather, 
some form of degenerate perturbation theory is used to 
bring the Hamiltonian into block diagonal form in dou- 
ble occupan cy space. More precisel y, one can perform 



a similaritv (|Eskes et all Il994i : iMacDonald et all Il988f ) 



transformation to bring the Hubbard model into block- 
diagonal form, in which the energies of a given block are 
approximately (in the limit of [/ 3> i) some number times 
U . Since the Hamiltonian is now block diagonal, it makes 
sense to project to the lowest energy block. The resultant 
Hamiltonian is, at leading order, the simpler model 



tj 



u ^ 



(17) 



known as the t — 3 model. Here hi — X^jo- Kt^zo-^jct where 
J is summed over the nearest neighbors of i and — 
— Vj, = 1. At half- filling, the first term vanishes as well 
as the three-site hopping process from the second term. 
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The result is a pure spin Hamiltonian 



(18) 



where each spin operator is a product of fermion bi- 
linears, S = 5*^* with 



Such a spin model is clearly invariant under the trans- 
formation 5* where h is an SU(2) matrix with 
unit norm. By contrast, the Hubbard model lacks this 
local SU(2) symmetry. It has only a global SU(2) as a 
result of the hopping term which is present even at half- 
filling. While it is certainly possible for a low-energy the- 
ory to possess different symmetries than the UV starting 
point, in this case the local SU(2) symmetry is spuri- 
ous as it appears entirely because we have ignored dou- 
ble oc cupancy. This discrepancy has be e n poi nted out 



before! Affleck et all 119881: iDagotto et all 119881 ) but the 



low-energy replacement for the Heisenberg model that 
preserves the global SU(2) symmetry of the Hubbard 
model was not offered. 

Strictly speaking, however, the block diagonalization 
procedure does permit double occupancy of bare elec- 
trons. This fact is hidden by the block diagonalization 
procedure itself because the resultant eigenstates, upon 
block diagonalization, are complicated linear combina- 
tions of the original electronic states. Consequently, the 
operators appearing in the t — J model are not the bare 
electrons in the Hubbard model. There are t/U correc- 
tions that fundamentally change the physics and account 
for i/ri/j > 1. As a result, a hole in the t — J model is 
not equivalent to a hole in the Hubbard model. Conse- 
quently, to extract the information that is in the Hub- 
bard model with the simpler t — J model, the operators 
must be re-transformed to the original basis which cer- 
tainly does not respect the no-double occupancy condi- 
tion. This procedure is cumbersome. For example, if one 
were to substitute the transformed electron operators in 
terms of the bare electrons of the Hubbard model into the 
Eq. ((T7)) . one would obtain Eq. (HUl) (see Appendix) as 
the true low-energy theory of the Hubbard model. While 
the first two terms in Eq. (|40p constitute the t — J model 
including the three-site hopping process, the remainder 
of the terms do not preserve the number of doubly occu- 
pied sites. As expected, the matrix elements that connect 
sectors which differ by a single doubly occupied site are 
reduced from the bare hopping t to /U . All such terms 
arise from the fact that the transformed and bare electron 
operators differ at finite U . As the terms involving double 
occupancy have the same amplitude as do the pure spin 
terms, there is no justification for dropping them. Hence, 
although t/U is small, (of order 1/10), it is sufficiently 
big to have important consequences which influence ex- 
perimental observables. The payoff is that hidden in the 



t/U corrections are emergent low-energy dynamics that 
are qualitatively different from that of the hard projected 
t-J m odel. In essence, there is a problem in the order of 
limits ljChov and Phillipsl . 120051 ) . The limits U ^ oo and 
i — > oo do not commute. That is, the physics is sensitive 
to a finite length scale for encountering a doubly occupied 
site. A comparison is helpful here between the predic- 
tions of the t-J and Hubbard models. The only limit in 
which the t-J model can be solved exactly is at the special 
value of the coupling J — t i n one spatial dimension i n 
which supersymmetry obtains! Bares and Blatter . IQQOl) . 
At this point, we can compare directly with the Hubbard 
model to see if anything is lost by ignoring double occu- 
pancy. The most striking differences are those for the 
exponents governing the asymptotic fall-off of the den- 
sity correlations, ac, and momentum distribution func- 
tions (9) in the t — J (with t = J, the supersymmet- 
ric point) and Hubbard models in d = 1. Here these 
quan t ities can be obtained e xactlvjFrahm and Korepin , 
19901 : iKawakami and Yangl . Il990a||d) for both models 
using Bethe ansatz. In the d=l Hubbard r riodel , 
the exponent 9 approaches ! Frahm and Korepin . IQQOt 
Kawakami and Yang , IQQObI) 1 /8 asymptotically as f/ — > 
oo fo r any filling. By contrast in the t — J 
model! Kawakami and Yan i, Il990ah . 9 is strongly de- 



pendent on doping with a value of 1/8 at half-filling 
and vanishing to zero as n decre ases. More surprising , 
the ex ponent ac remains pinned !jKawakami and Yang . 
1990allbi ) at 2 for the t/ -)> oo limit of the Hub- 
bard mode l at any filling. In fact, a t any value of 
U, ac = 2! Kawakami and Yanj . Il990al lbl) in the dilute 
regime of the Hubbard model in d =1. In the t — J 
model! Kawakami and Yand . ll990al ) (t = J), ac starts at 
2 at n = 1 and approaches a value of A aX n — 0. Note, 
Q!c = 4 is the non-interacting value. That is, in d = 1 in 
the dilute regime, density correlations decay as r^^ in the 
[/ — > oo Hubbard model and as r"'' in the t — J (t=J). 
This discrepancy is a clear indicator that relevant low- 
energy physics is lost if double occupancy of bare elec- 
trons is projected out in the parameter range considered 
here. Supposedly, this is captured by the t/U corrections 
to the electron operators in Eq. !(39|). While it might 
not be surprising that t/U corrections are important for 
U/t = 2 as in the supersymmetric t-J model, similar dis- 
crepancies have been noted even for U/t ^ 10 in a direct 
comparison between the density o f states of the t-J a nd 



Hubbard models. The conclusion !|Leung et al\ . 119921 ) of 



such studies is that even at strong coupling values of 
U/t, the t-J and Hubbard models yield similar density of 
states for the lower Hubbard band only at dilute fillings 
starting around n = 1/2, considerably far away from the 
n — 1 point relevant to the cuprates. Perhaps the limi- 
tations of the t — J model in the context of the cuprates 
are best sum marized by one o f its chief practioneers, P. 
W. Anderson !jAndersonl, 1997 ): "From this point of view, 
it is nonsense to consider J values which are larger than 
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some small limiting value. All the theoretically excit- 
ing possibilities proposed by various advocates of the t-J 
model appear in this unphysical regime and relate to no 
real physical system. The spate of gauge theories, phase- 
separation theories, etc. which have their existence in 
this large-J region are therefore almost devoid of physi- 
cality." 



B. Exact Theory: Charge 2e Boson 



The remedy is to perform the Wilsoninan (| Wilson 



19711 ) procedure exactly for the simplest model of a doped 
Mott insulator that captures dynamical spectral weight 
transfer, namely the Hubbard model. The result will be 
a theory that contains all the degrees of freedom needed 
to capture the low-energy spectrum of a doped Mott in- 
sulator. Indeed it is dynamical spectral weight transfer 
that makes the construction of a low-energy theory of a 
doped Mott insulator non-trivial. In the exact theory, 
such effects are mediated by a collective charge 2e bo- 
son. The procedure is as follows. Consider hole-doping a 
Mott insulator. The high energy scale is now the up- 
per Hubbard band and hence must be integrated out 
to acquire the true low-energy theory. We introduce an 
elemental field which represents the creation of excita- 
tions in the upper Hubbard band. I n this regard, the 



19861) procedure 



slave-boson (jKotliar and Ruckenstein , 
has been used. The key idea here is to construct a model 
which has the same matrix elements as does the Hub- 
bard model in which the interaction t erm is replaced by 
[/dldl Kothar and Ruckenstein . 19861 ). where di is purely 
bosonic. One cannot do this, however, withou t simul- 
taneously introducing (|Kotliar and RuckensteiJ . Il986l ) 
three other bosonic fields and three other Lagrange mul- 
tipliers. The low-energy theory is obtained by integrating 
over the di field as it has a mass U. Since the di field re- 
ally does represent the creation of double occupancy, in- 
tegrating it out also gets rid of double occupancy. Hence, 
this procedure does not retain dynamical spectral weight 
transfer, the very feature we are trying to isolate by the 
Wilsonian procedure. Simply, this procedure does too 
much. Although auxiliary fields are present that repre- 
sent strictly singly occupied physics, there is no doube 
occupancy left when di is integrated out. We outline 
a procedure that does precisely what is demanded on a 
Wilsonian account but no more. 

The solution is to expand the Hilbert space so that the 
field associated with the UHB, and hence with mass U, 
represents the creation of double occupancy only when 
a constraint is imposed. Unconstrained, it mediates 
whatever physics transpires in the UHB. The true low- 
energy theory will then correspond to integration over 
the new field without imposing the constraint. As a 
consequence, the low-energy theory will explicitly per- 
mit double occupancy and hence have the essential in- 



gredients to describe dynamical spectral weight trans- 
fer. In so doing, the theory will contain the constraint 
parameter which permitted an identification of the new 
field with the creation of double occupancy in the first 
place. This pr ocedure is analogous t o that used by 
Bohm and Pines (|Bohm and PinesI , Il953l ) in their deriva- 
tion of pl asmons. In the i r cla ssic derivation, Bohm 



and Pines (|Bohm and PinesI . 119531 ) expanded the Hilbert 
space to include a collection of low-energy oscillators 
but a constraint was imposed such that when it was 
solved, the original electron gas Hamiltonian was recov- 
ered. When the constraint was relaxed, the plasmon 
emerged as a new propagating degree of freedom. Our 
construction mirrors theirs in essence. 

To this end, we extend the Hilbert space of the Hub- 
bard model to include on each site a new fermionic oscil- 
lator, Di, which will represent the physics of the UHB. 
Through a constraint, will represent the creation of 
double occupancy. Imposing such a constraint requires 
a trick because double occupancy transforms as a boson 
as it involves the product of two fermionic operators. At 
the same time, there can only be one double occupancy 
per site. Hence, double occupancy has both fermionic 
and bosonic character. Dealing with this dual character 
requires a trick from supersymmetry. Let us recall how 
supersymmetry was first introduced into string theory. 
In its earliest form, string theory was a theory entirely of 
bosonic modes represented by some field (a) , where a 
is the distance along the string and /i are the cartesian 
coordinates. Ultimately to make this theory applicable 
to anything real, such as QCD, fermio ns had to be in- 
cluded. To solve this problem, Ramond (| Ramondl . Il97l[ ) 
promoted the Clifford matrices to the level of a field and 
defined a super-field 



(20) 



where is a non-commuting complex Grassmann param- 
eter and 7/i(a) represent the Clifford matrices. Since 
6*^ — 0, the second term in Xj^{a,9) transforms as a bo- 
son. Hence, fermions can be included in a bosonic theory 
simply by multiplying a fermion by a Grassmann. In 
our problem, just the opposite procedure will be used. 
We can turn a boson, namely double occupancy, into a 
fermion by multiplication with a Grassmann. The con- 
straint for the elemental field Di will turn out to be 
roughly S{Di — Oci^Cii). However, our theory is not 
supersymmetric as Di cannot properly be regarded as 
a superfield because it really does not mix bosons and 
fermions as does Xfj_{a,9). The introduction of the su- 
persymmetric coordinate 9 is simply a trick to integrate 
out the high-energy physics. Nowhere in the UV or IR 
limits of the theory will 9 appear. This is taken care 
of by writing the Lagrangian explicitly as an integration 
over Grassmans. To construct the Lagrangian in this ex- 
tended space, one has to include the standard dynamical 
terms for the elemental fields, Ci^ and Di, as well as the 
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hopping terms such that when the constraint is solved, 
one recovers the Hubbard model. The Lagrangian that 
fits the bill is 



- {D]0Ci^aVcrCj-cr + h.C 



i 



(21) 



where Cij^a = 



aij^a = d6{\ — ni^^cr){^ — nj^^a) and dPO 
represents a complex Grassniann integration. In this La- 
grangian, the first two terms represent the dynamics of 
electrons in the lower Hubbard band and the heavy field, 
respectively, the third electron hopping in the LHB, the 
fourth decay of the heavy field into an electron singlet 
state on neighbouring sites and the fourth an the creation 
of heavy field particle-hole pair accompanied by the cre- 
ation of an electron-hole pair on neighbouring sites. The 
constraint Hamiltonian i/con is taken to be 



0c 



) + h.C. 



(22) 



and arises simply by exponentiating a (5— function. The 
constant s has been inserted to carry the units of energy. 
The precise value of s will be determined by compar- 
ing the low-energy transformed electron with that in Eq. 
([39)l . This Lagrangian was constructed so that if we solve 
the constraint, that is, integrate over if and then Di, we 
obtain exactly / d'^OOOLnuhh = -^Hubb, the Lagrangian of 
the Hubbard model. Hence, up to a factor of unity, our 
starting Lagrangian is equivalent to the Hubbard model. 

The advantage of our Lagrangian is that it permits a 
clean identification of the f7-scale physics without equat- 
ing it with double occupancy. As the theory is Gaussian 
in the massive field, it can be integrated out exactly. To 
accomplish this, we define the matrix 



(23) 



and bi = I]j &u = Y,ja9ijCj,crVcrCi.a- We now com- 
plete the square and integrate over Di exactly in the 
partition fu nction. At zero fr equency the low-energy 
Hamiltonian (|Chov et all l2008al ) is 



where 



U 



2J 



'E "^^J'tc.a + f E v\{M-%jb, + h.C. (24) 



which constitutes the true (IR) limit as the high-energy 
scale has been removed. This Hamiltonian appears com- 
plicated but can be handled simply as will be seen. In 
essence, the first term in iJint, which contains the spin- 
spin interaction and three-site hopping terms, is irrele- 
vant relative to the terms linear in b so far as the charge 
excitations are concerned. It is this reduction that makes 
it possible to isolate the propagating degrees of freedom 
in a doped Mott insulator. 

To fix the energy scale s, we determine how the electron 
operator transforms in the exact theory. As is standard, 
we add a source term to the starting Lagrangian which 
generates the canonical electron operator when the con- 
straint is solved. For hole-doping, the appropriate trans- 
formation that yields the canonical electron operator in 
the UV is 

-n)cL + v„D]ec,_^.. 



i.a 



99(1 



However, in the IR in which we only integrate over the 
heavy degree of freedom, Di, the electron creation oper- 



ator becomes 

J 



(1 - n 



r)c. 



,t 

i.a 



^ t 



U 



(25) 



to linear order in t/U. This equation bares close resem- 
blance to the transformed electron operator in Eq. p9p . 
as it should. In fact, the first two terms are identical. 
The last term in Eq. ([39]) is associated with double oc- 
cupation. In Eq. ([25l) . this role is played by (pi. De- 
manding that Eqs. (|39p and ([25|l agree requires that 
s = t, thereby eliminating any ambiguity associated with 
the constraint field. Consequently, the complicated in- 
teractions appearing in Eq. (l40l) as a result of the in- 
equivalence between the tranformed and bare fermions 
are replaced by a single charge 2e bosonic field (pi which 
generates dynamical spectral weight transfer across the 
Mott gap. While it is traditional in solid state systems 
in which both bosons and fermions appear to integrate 
out the bosons, that would be incorrect here. Both the 
bosons and fermions are light degrees of freedom that 
mediate low-energy physics. 

While electron number conservation is broken in the 
IR, we find by inspection of Eq. (p4)) that a conserved 
low-energy charge does exist, given by 



+ 2 



(26) 



2, J 



Physically, Q should equal the total electron fill- 
ing, implying immediately implying that the weight 
of the low-energy fermionic part must be less than 
the conserv ed charge. In fact, E q. ((26)) gives a 
prescription! Chakraborty et al. . 20091 ) for a, namely the 
bosonic charge, if we interpret Q as 1 — a; and the 
fermionic quasiparticle density as 1 — x', thereby cor- 
roborating the quasiparticle picture in Fig. (|5Jd). 
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1. Half-Filling: Bound doublon/holon pairs 

As it is the dynamical part of the spectral weight that 
is intimately connected with the pseudogap, it stands to 
reason that the ultimate role of the charge 2e boson is 
to mediate bound charge excitations. Physically, it is 
reasonable that the charge 2e boson can only influence 
the dynamics at low energies through the formation of 
bound states because it does not have a Fock space of 
its own. That is, once the heavy field is integrated out, 
the Hilbert space is that of the Hubbard model. Imme- 
diate evidence that this state of affairs, namley bound- 
state formation, obtains arises from the equivalent theory 
at half-filling. At half-filling, both the upper and lower 
Hubbard bands reside at high energy and hence must be 
integrated out. This can be accomplished by introduc- 
ing an additional fermionic field that represents the cre- 
ation of excitations in the LHB. Its associated Lagrange 
multiplier will be (pj. Unlike the theory at finite dop- 
ing, t he theorvdChoy et al. . 2008bt Leigh and Phillipsl . 



20081 : 

the pesky 



Phillips et al 



)y et al 
120091 ) 



at half-filling will not have 
matrix, A4 and as a result, no 
The exact low- e nergy 



fermionic 

bare fields will have dynamic s ^„ 

Lagrangian (|Chov aLl . l2008bHPhillips et aLl . l2009[) that 
arises from this procedure, 

2 



rhf _ r, 



S 

U ' ~ U 



(27) 



•,\l'^\<./2+p,ui/2+ui' ,-\C-\i/2-p,u>/2~ui' ,i 
+ h.C, (28) 

contains two bosonic fields with charge 2e ((^^) and 
— 2e {(p). These bosonic modes are collective degrees 
of freedom, not made out of the elemental excitations. 
They represent dynamical mixing with scale physics, 
namely, the contribution of double holes (— 2e) and dou- 
ble occupancy (2e) to any state of the system. Here s 
and s are constants with units of energy, all operators in 
Eq. ([27| have the same site index, repeated indices are 
summed both over the site index and frequency, w, c|g. 
creates a fermion on site i with spin cr. 



- ^ e^p'' Ck/2+p,tCk/2-p,i, 



(29) 



-P 



and the dispersion is given by Ep 

cos(A:^a/2) cos(p^a), where k and p are the 
center of mass and relative momenta of the fermion pair. 
The coefficients 



■ 2a; 



i'^' (-) = + + v/T^Mf/ (30) 

play a special role in this theory as they account for the 
turn-on of the spectral weight. At the level of a La- 



grangian, the vanishing of the coefficient of a quadratic 
term defines the dispersion of the associated particle. All 
the terms which are naively quadratic, Eq. (f27|) . possess 
constant coefficients and hence we reach the conclusion 
that there are no bare propagating bosons or electrons. 
In fact, it is the vanishing of the pesky calM matrix in 
the half-filled system that leads to this state of affairs. 
This implies that there is no spectral weight of any kind. 
However, a Mott insulator has spectral weight as depicted 
in Fig. (j4]). Consider the second line of the Lagrangian, 
Eq. (j28|) . Appearing here are two interaction terms, 
which describe composite excitations, whose coefficients 
can vanish. These operators might then be thought of 
as the kinetic terms for composite excitations mediated 
by the charge ±2e bosonic fields (loosely speaking, we 
might think of this as occurring because of the formation 
of bound states). Such an interpretation is warranted be- 
cause the spin-spin interaction and all higher-order oper- 
ators contained in the |&[^ term are at least proportional 
to a'* and hence are all sub-dominant to the composite 
interaction terms. Here a is the lattice constant. Con- 
sequently, at the level of the Lagrangian, the turn-on 
of the spectral weight is governed by the vanishing of 
the coefficients of the coupled boson-fermion terms. Fig. 
([M)) shows explicitly that the vanishing of 7 and 7 leads 
to spectral weight which is strongly peaked at two dis- 
tinct energies, ±C//2. Each state in momentum space 
has spectral weight at these two energies. The width of 
the bands is 8t. The particles which give rise to the turn- 
on of the spectral weight are composite excitations or the 
bound states of the bosonic and fermionic degrees of free- 
dom determined by the interaction terms ip'^cc and (pec. 
In the terms of the variables appearing in the Hubbard 
model, the composite excitations correspond to bound 
states of double o ccupancy and holes as has been pos- 
tulated previouslv (|Castellani et to be the ulti- 
mate source of the gap in a M ott insulator. Interestin gly, 
the slaved-particle approach ( Castellani et al. . 1992f) on 



the large N limit of the infinite- 1/ Hubbard model also 
finds that the Mott gap originates from a gap in the spec- 
trum of an auxiliary boson. In so far as they generate 
the spectral weight, the interaction terms can be thought 
of as the kinetic terms in the low-energy action. The 
gap (Mott gap) in the spectrum for the composite ex- 
citations obtains for U > 8i as each band is centered 
at ±/7/2 with a width of 8t. Fig. (fT4|) demonstrates 
that the transition to the Mott insulating state found 
here proceeds by a discontinuous vanishing of the spec- 
tral weight at the chemical potential to zero but a con- 
tinuous evol ution of the Mott gap as is seen in numerical 
calculations (IPark et all 



2008f) in finite-dimens ional lat- 
o ljCeorges et al. , 1996[ ) solution. 



tices but not in the d = 

In terms of the bare electrons, the overlap with the 
composite excitations determines the Mott gap. To de- 
termine the overlap, it is tempting to complete the square 
on the (^^cc term bringing it into a quadratic form, ^t^* 
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with ^ = A(p + Bcc. This would lead to composite exci- 
tations having charge 2e, a vanishing of the overlap and 
hence no electron spectral density of any kind. However, 
the actual excitations that underlie the operator ip'^cc 
correspond to a linear combination of charge e objects, 
and V3^c. In terms of the UV variables, the latter 
can be thought of as a doubly occupied site bound to 
a hole. To support this claim, we re-do the procedure 
quoted ab ove for generating t he electron operator . At 
half-filling (jChojLelo/J, I20n8bl : IPhillips et a/1 120091 ). the 
exact representation of the electron creation 



u 



(31) 



is indeed a sum of two composite excitations, the first 
having to do with spin fiuctuations (6^c) and the other 
with high-energy physics, if^c and (pc, that is, excitations 
in the UHB and LHB, respectively. We can think of the 
overlap 



0= |(ct|5t)(gt|^t)|2p^ 



(32) 



in terms of the physical process of passing an electron 
through a Mott insulator. The overlap will involve that 
between the bare electron with the low-energy excitations 
of Eq. (|31]) . (c|c), and the overlap with the propagating 
degrees of freedom, {c\^) with P^, the propagator for 
the composite excitations. Because of the dependence 
on the bosonic fields in Eq. ([3T|) . O retains destructive 
interference between states above and below the chemi- 
cal potential. Such destructive interference between ex- 
citations across the chemical potential l eads to a vanish- 



ing o f the spectral weight at low energies (jMeinders et al. 
19931) . Consequently, the turn-on of the electron spectral 



weight cannot be viewed simply as a sum of the spec- 
tral weight for the composite excitations. As a result of 
the destructive interference, the gap in the electron spec- 
trum will always exceed that for the composite excita- 
tions. Hence, establishing (Fig. (HU) that the composite 
excitations display a gap is a sufficient condition for the 
existence of a charge gap in the electron spectrum. A sim- 
ple calculation. Fig. ([M]). of the electron spectral func- 
tion a,t U — St confirms this basic principle that a gap in 
the propagating degrees of freedom guarantees that the 
electron spectrum is gapped. Further, Fig. p4)) confirms 
that the electron spectral function involves interference 
across the Mott scale. Consequently, although the com- 
posite excitations are sharp, corresponding to poles in 
a propagator as in Eq. (1^51) . the electrons are not. In 
fact, as Eq. (j3ip lays plain, an electron is in a linear 
superposition of excitations both in the lower and upper 
Hubbard bands. Consequently, in terms of the original 
electron degrees of freedom, the transition to the Mott 
gap will involve spectral weight transfer at energies on 
the Mott scale as illustrated in Fig. (jS]). Alternatively, 
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FIG. 14 a) Diamond-shaped surface in momentum space 
where the particle dispersion changes sign, b) Turn-on of the 
spectral weight in the upper and lower Hubbard bands for the 
composite excitations as a function of energy and momentum. 
In the UHB, the spectral density is determined to 7p while for 
the LHB it is governed by 7p. The corresponding operators 
which describe the turn-on of the spectral weight are the com- 
posite excitations ip^ cc (UHB) and (pec (LHB). The electron 
spectral density is determined by an overlap (see Eq. (I32p ) 
with these propagating collective modes, c) Spectral function 
for the electrons aX U = 8t. Clearly shown is the gap in the 
spectrum and non-zero spectral weight at all momenta in the 
first Brillouin zone. 



if the experimental probe were the composite or bound 
states, spectral weight transfer would be absent because 
the composite excitations represent the orthogonal prop- 
agating modes of a half-filled band. 

The composite excitations found here described by the 
vanishing of 7p and 7p are the propagating degrees of 
freedom in a Mott insulator. They are orthogonal in the 
sense that they never lead to a turn-on of the spectral 
weight for the composite excitations in the same energy 
range. This analysis demonstrates that the spin-spin in- 
teraction, contained in the |6p term, plays a spectator 
role in the generation of the Mott gap. Nonetheless, 
there is a natural candidate for the antiferromagnetic or- 
der, namely Bij — {gijLp\ci^^Cj_\) . The vacuum expec- 
tation value of this quantity is clearly non-zero as it is 
easily obtained from a functional derivative of the parti- 
tion function with respect to 7p. Such an antiferromag- 
net has no continuity with that of weak-coupling theory. 
Hence, both the Mott gap and subsequent antiferromag- 
netic order emerge from composite excitations that have 
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no counterpart in the original UV Lagrangian but only 
become apparent in a proper low-energy theory in which 
the high-energy degrees of freedom are explicitly inte- 
grated out. In the doped state, a similar state of affairs 
obtains. 



2. Experimental Consequences 



As Eq. ([25|) lays plain, the fermionic operator at 
low energies is a linear superposition of two excitations. 
The first is simply the standard excitation in the LHB, 
(1 — ni^-a-)cl ^ {ni^-aCia m thc UHB for electron dop- 
ing) with a renormalization from spin fluctuations (sec- 
ond term). The second is a new charge e excitation, 
Ci-aM~j^(p'j. In the lowest order in t/U, our theory 

predicts that the new excitation corresponds to Ci^-afl, 
that is, a hole bound to the charge 2e boson. This extra 
charge e state mediates dynamical (hopping-dependent) 
spectral weight transfer across the Mott gap. A saddle- 
point analysis will select a particular solution in which 
ifi is non-zero. This will not be consistent with the gen- 
eral structure of Eq. ([25]) in which part of the electronic 
states are not fixed by ipi. Similarly, mean field theory in 
which ipi is assumed to condense, thereby thwarting the 
possibility that new excitations form, is also inadequate. 
The procedure outlined in the second Appendix circum- 
vents these problems and preserves the integrity of Eq. 
(1251). 

The resultant spectral function U — IQt is displayed in 
Figs. ([T5|) and (flB)) . First, a low-energy kink is present in 
the electron dispersion for a wide range of doping. In fact, 
more than one kink exists as is evident from the enlarged 
region. Fig. (fTSk): 1) one at roughly Q.2t « lOOmeV 
and the other at 0.5t « 250meV. To pin-point the origin 
of these kinks, we treated the mass term of the boson 
as a variable parameter and verified that the low-energy 
kink is determined by the bare mass. In the low-energy 
theory, the bare mass of the boson is t'^/U, independent 
of doping. Both the doping dependence and energy of 
this k ink are consistent with experiment (iLanzara et al 



20011 ) . Wh ile phonons(|Lanzara et ali. l200lh an d spin 



fluctuations (|Macridin et all l2007t IZemliic et all l2008l ) 



have been invoked to explain the low-energy kink, the 
hidden charge 2e boson offers a natural explanation 
within the strong correlation physics of the Mott state. 

At sufhciently high doping (see Figs. (fTHk) and (fTBb')')- 
the high-energy kink disappears. Experimentally the ori- 
gin o f the high-energy k ink is hotly debated. In fact, 
some dnosov et al. . doubt its intrinsic importance, 

attributing it to an extraneous matrix element ef fect. In 



the initial experiments by Graf, et al. (|Graf et a/.l . l2007f ). 



the high-energy kink is accompanied b y a splitting of th e 



electron dispersion into two branches (jOraf et all 120071 ). 
The two branches were interpreted as evidence for spin- 
charge separation. As evident from Fig. ((T5|) . the high- 
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FIG. 15 (a) Spectral function for filling n = 0.9 along the 
nodal direction. The intensity is indicated by the color 
scheme, (b) Location of the low and high energy kinks as 
indicated by the change in the slope of the electron disper- 
sion, (c) The energy bands that give rise to the bifurcation 
of the electron dispersion. 



energy kink is proceeded by a bifurcation of the elec- 
tron dispersion below the chemical potential into two 
branches. 

The energy difference between the two branches 
achieves a maximum at (0, 0) as is seen experimentally. 
A computation of the spectral function at U = 20t and 
n = 0.9 reveals that the dispersion as well the bifurcation 
still persist. Further, the magnitude of the splitting does 
not change, indicating that the energy scale for the bi- 
furcation and the maximum energy splitting are set by t 
and not U. The origin of the two branches is captured in 
Fig. ((T5b ). The two branches below the chemical poten- 
tial correspond to the standard band in the LHB (open 
square in Fig. (jl5b ) on which ip vanishes and a branch 
on which (f ^ (open circles in Fig. ([T5b). The two 
branches indicate that there are two local maxima in the 
integrand in Eq. (|Tf|. One of the maxima, tp ~ 0, arises 
from the extremum of G{k,Lu,ip) whereas the other, the 
effective free energy (exponent in Eq. (j47l) ) is minimized 
{if ^ 0). Above the chemical potential only one branch 
survives. The split electron dispersion below the chem- 
ical potential is consistent with the composite nature of 
the electron operator dictated by Eq. (^5]) . At low ener- 
gies, the low-energy fermions are linear superpositions of 
two states, one the standard band in the LHB described 
by excitations of the form, cl^{l ~ Uig) and the other 
a composite excitation consisting of a bound hole and 
the charge 2e boson, Cisip\. The former contributes to 
the static part of the spectral weight transfer (2x) while 
the new charge e excitation gives rise to the dynamical 
contribution to the spectral weight transfer. Because the 
new charge e state is strongly dependent on the hopping 
it should disperse as is evident from Fig. ([T5)) and also 
confirmed experimentally. The operator that describes 
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FIG. 16 Spectral function for two different fillings (a) n = 0.8 
and (b) n = 0.4 along the nodal direction. The absence of 
a splitting in the electron dispersion at n — 0.4 indicates 
the bifurcation ceases beyond a critical doping. The spec- 
tral functions for two different values of the on-site repulsion, 
{c)U = lot and {d)U = 20t for n = 0.9 reveals that the high- 
energy kink and the splitting of the electron dispersion have 
at best a weak dependence on U. This indicates that this 
physics is set by the energy scale t rather than U. 



this excitation is given by Eq. (|45|) and it constitutes the 
propagating degree of freedom in a doped Mott insula- 
tor as it generates a pole in the composite particle Green 
function (Eq. (|48| ). In the composite basis, the electron 
are not sharply defined. 

Because the charge 2e boson is a local non-propagating 
degree of freedom, the formation of the composite exci- 
tation, Ciff(p'' leads to a pseudogap at the chemical po- 
tential. The spectral functions. Fig. aX n — 0.9 
and n = 0.8 reveal an absence of spectral weight at the 
chemical potential. In the strongly overdoped regime, 
spectral weight emerges at the chemical potential and the 
pseudogap vanishes. The formation of a gap in a single- 
band system, requires a vanishing of single-particle Green 
function. This is mediated by a divergence of the elec- 
tron self-energy along a connected surface in momentum 
space. Computed in Fig. (fT7| is the imaginary part of 
the self energy at different temperatures. At low temper- 
ature (T < f^/U), the imaginary part of the self-energy 
at the non-interacting Fermi surface develops a peak at 
u! = 0. At T — 0, the peak leads to a divergence. This is 
consistent with the opening of a ps e udoga p. As we have 
pointed out earlier! Stanescu et al . 2007t ). a pseudogap 
is properly identified by a zero surface (the Luttinger 
surface) of the single-particle Green function. This zero 






> 1 ' 




- T=0.05t 






- T=O.U 






T=0.2t 






- T=0.4t 





2 
oo/t 



FIG. 17 The imaginary part of the self energy as the function 
of temperature for n = 0.7. A peak is developed at = 
at low temperature which is the signature of the opening of 
the pseudogap. The density of states explicitly showing the 
pseudogap is shown in adjacent figure. 



surface is expected to preserve the Luttinger volume if 
the pseudogap lacks particle-hole symmetry as shown in 
the second of the figures in Fig. PT| . 

Using this approach, we are also able to address the 
origin of the mid-infrared band (MIB) in the optical con- 
ductivity. From the Kubo formula. 



<yxxiuj) = 27re 



(Pk J duj'{2tiiink^f 



A{uj +uj' ,k)A{uj\k) 
(33) 

we computed the optical conductivity. Here /(w) is the 
Fermi distribution and A(w, k) is the electron spectral 
function. At the level of theory constructed here, the 
vertex corrections are all due to the interactions with 
the bosonic degrees of freedom. Since the boson acquires 
dynamics only through electron motion and the leading 
such term is 0{t^ /U^), the treatment here should suffice 
to provide the leading behaviour of the optical conduc- 
tivity. Further, to isolate the mid-infrared, we subtracted 
the Drude weight at the origin. Apparent in the optical 
conductivity shown in Fig. ([T8| is a peak at ui/t « .5t. 
This constitutes the mid-infrared band. From the inset, 
we see that the frequency, Wmax at which the MIB ob- 
tains is an increasing function of the electron function 
filling, whereas the integrated weight, N^si^), is a de- 
creasing function as the filling increases, both of which 
are in agreement with experiment. We set the integration 
cutoff to = 2< = 1/m*. However, it does not vanish 
at half-filling. Experimentally, A^eff also does not vanish 
when extrapolated to half-filling. The persistence of the 
the MIB at half-filling suggests that the mechanism that 
causes the MIB is apparent even in the Mott state. We 
determined what sets the scale for the MIB by studying 
its evolution as a function oi U. As is clear from Fig- 
ure (|18p . Wmax is set essentially by the hopping matrix 
element t and depends only weakly on J. The physical 
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FIG. 18 Optical conductivity as a function of electron filling, 
n, with Drude weight at origin subtracted. The peak in the 
optical conductivity represents the mid-infrared band. Its ori- 
gin is mobile double occupancy in the lower-Hubbard band. 
The insets show that the energy at which the MIB acquires 
its maximum value, ojmax is an increasing function of electron 
filling. Conversely, the integrated weight of the MIB decreases 
as the filling increases. This decrease is compensated with an 
increased weight at high (upper-Hubbard band) energy scale. 



processes that determine this physics are determined by 
the coupled boson-Fermi terms in the low-energy theory. 
The (filci^Cii term has a coupling constant of t whereas 
the (flbi scales as t'^/U. Together, both terms give rise to 
a MIB band that scales as Wmax/^ = 0.8 — 2.21t/U (see 
inset of Fig. (dJ)). Since t/U ^ 0(.l) for the cuprates, 
the first term dominates and the MIB is determined pre- 
dominantly by the hopping matrix element t. Within the 
interpretation that ip represents a bound state between 
a doubly occupied site and a hole, second order pertur- 
bation theory with the flbi term mediates the trans- 
port of hole over two sites mediated by an intervening 
doubly occupied site (see Fig. (1) of Ref. ljChov et all 
2008al )). It is the resonance between these two states 
that results in the mid-IR band. Interestingly, this res- 



onance persists even at half-filling and hence the non- 
vanishing of iVoff at half-filling is not evidence that the 
cuprates are n ot doped Mott insulat ors as has been re- 
cently claimed (|Comanac et aU 120081 ). Rather the quan- 
tum fluctuations that are present even in the half-filled 
system still persist at finite doping. In the half-filled sys- 
tem, bound states form which provide a gap to single 
particle excitations. In the doped system, such excita- 
tions are only partially gapped (in momentum space as 
evidenced by the V-shaped gap in Fig. dTTt ). giving rise 
to a pseudogap. 

One of our key contentions is that the correct theory 
of the pseudogap should at high temperatures explain 
T— linear resistivity in the strange-metal regime. Indeed 
this theory has the ingredients to do this. The mecha- 
nism is simple. In the pseudogap regime, the charge 2e 
boson is bound. This gives rise to an explicit violation 
of the band-insulator sum rule that L = Uh- Once the 
T* line is crossed, Lp unbinds. Hence, above the T* line. 
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FIG. 19 Proposed phase diagram for the binding of the holes 
and bosons that result in the formation of the pseudogap 
phase. Once the binding energy vanishes, the energy to ex- 
cite a boson vanishes. In the critical regime, the dominant 
scattering mechanism is still due to the interaction with the 
boson. T- linear resistivity results anytime T > cot, where uub 
is the energy to excite a boson. To the right of the quantum 
critical regime (QCP), the boson is irrelevant and scattering 
is dominated by electron-electron interactions indicative of a 
Fermi liquid. The QCP signifies the end of the binding of 
fermi and bosonic degrees of freedom that result in the pseu- 
dogap phase. 



the temperature exceeds the energy to create the charge 
2e boson. However, it still scatters off the electrons. The 
problem is now a trivial one of electron-boson scattering 
above the temperature to create the boson. This is an 
old problem and the result is well known. The resistiv- 
ity scales as a linear function of temperature. Hence, as 
depicted in Fig. p9|) . the strange metal emerges as the 
unbound phase of the charge 2e boson in which critical 
fluctuations are the scattering mechanism. Above the T* 
line, L = Uh- Hence, this mechanism makes a clear ex- 
perimental prediction that the T* line is the boundary 
below which L > 2x obtains. A repetiti on of the experi- 
ments of Chen, et al. ljChen et al ]. ll99ll) as a function of 
temperature would directly falsify this claim. Further, 
the mechanism for the strange metal regime is consistent 
with the implication of the scaling analysis leading to Eq. 
(jlip . Namely, T— linear resistivity requires an additional 
energy scale absent from a single-parameter scaling anal- 
ysis. In the exact low-energy theory, a charge 2e boson 
emerges as a new degree of freedom. While it is bound 
in the pseudogap regime, its unbinding beyond a critical 
temperature or doping provides the added degree of free- 
dom to generate the anomalous temperature dependence 
for the resistivity. A further experimental prediction of 
this work then is that the strange metal regime should be 
populated with charge 2e excitations, without the usual 
diamagnetic signal. Shot noise measurements are ideally 
suited for testing this prediction. 
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V. OUTLOOK AND PREDICTIONS 

The central claim in this Colloquium is that composite 
excitations, doublon-holon pairs, emerge as the propa- 
gating degrees of freedom in the normal state of a doped 
Mott insulator. In light of other strongly coupled prob- 
lems such as QCD, this state of affairs is not surprising. 
This physics arises because although double occupancy 
mixes into the ground state, the empty sites which are 
dynamically generated are not free to move around. Con- 
sequently, the low-energy physics in a doped Mott system 
is determined by the effective doping level x' = a; -t- a, 
where a reflects the dynamical hole count. Such physics 
is mediated through the charge 2e boson. At half-filling, 
the band structure of the doublon-holon pairs leads to 
the gapp ed sp e ctrum of a Mott i nsulator as antic i pated 
by MottdMottl. I1949I) and others (ICastellani et ^ . Il979t 
Kohnl . ll964l : lYokovama tt al\ . I2OO6I ). Because the excita- 
tions described by 7 and 7 never share a common energy 
where the spectral weight is non-zero, they can be viewed 
as the independent propagating degrees of freedom in a 
half-filled band. Both the metallic state at half-filling 
and the strange metal are mediated by the unbinding of 
the composite excitations. The simplest way of under- 
standing why the charge 2e boson must be bound at low 
energies, aside from the fact that it has no bare dynam- 
ics, is that once the high-energy sector is integrated out 
exactly, the Hilbert space shrinks back to the Fock space 
of the Hubbard model. The charge 2e boson acting in 
this Hilbert space can only mediate dynamics through 
binding with the elemental fields. 

A clear indicator that the theory pre- 
sented h ere is on the right track is the 
calculation! Chakrabortv and Philhpsl . I2OO9I) shown 
in Fig. ((20|) of the carrier density as revealed by the 
Hall number. Clearly shown are two components to 
the carrier density, one temperature independent the 
other a highly temperature dependent component which 
gives rise to activated behaviour in agreement with 
the phenomenological fit of the experimental data in 
Eq. dll). In fact, the gap, A (see inset of Fig. ([20])) is 
in excellent agreement with the experimental data on 
LSCO. Above T*, the gapped component vanishes. Such 
two-fluid behaviour suggests that in the effective doping 
level, a is the temperature-dependent component. Any 
experimental probe that couples to the low-energy 
excitations should be interpreted in terms of a;', not the 
bare hole number x. Several experimental predictions 
follow: 1) above a critical doping level, a should (as 
a result of the unbinding of the charge 2e boson) and 
L/uh should be doping dependent, alread y confirmed 



3) the superfluid density should exceed x and scale as as 
x-\-_a, already confirmed in YBa9Cu':iOfi-|--r l|CooDer et al. 



lionh (YBCO), 4) the inverse dielectric function should 



by re cent x-ray oxygen K-edge experiments (jPeets et al 



2009f ). 2) similarly, because the T* line corresponds to 
the unbinding of the charge 2e boson, angle-resolved 
photoemission experiments should observe a narrowing 
of line shapes as the temperautre is increased above T* , 



possess two-particle hole continual Chov et al . 2008al ) 



the second feature (starting at approximately 0.5i) 
reflecting the new bound state and 5) Fermi surface vol- 
umes, that is the total volume of the hole pockets minus 
that of the electron poc kets, extracted from q u antum 



oscillation experiments (jPoiron-Levraud et al. . 2007 ) 



should be compared with 2 (a; -I- a) not 2x. The latter is 
particularly germane beca use the Fermi surface v o lumes 



extracted experime ntally ()Doiron-Levraud et all 12007 



LeBoeuf et aU . l2007t l for YBCO are not consistent with 
any integer multiple of the physically doped holes. 

Theoretically, several questions remain. First, is the 
theory in the composite excitation picture natural in the 
sense that there are no relevant perturbations to normal- 
state physics. This would require a calculation of the 
/3-function which necessitates rewriting the pure electron 
and boson interactions in terms of the composite parti- 
cles and then a subsequent renormalization group analy- 
sis. At present it is unclear how to proceed along these 
lines. Second, can it be quantified at what doping level 
does the charge 2e boson decouple from the low-energy 
physics and a Fermi liquid description becomes valid. 
Third, what role do the composite excitations play in 
the superconducting state. Since these excitations arise 
from a mixing with the high-energy sector, they have a 
chance of accountng for the otherwise unexplained exper- 
imental observ ation that the onset of superconductivity 
in the cuprates ljMolegraaf et aLl . l2002n is accompanied by 
a depletion of spectral weight at high energies {U— scale 
physics) and a compensating increase at low energies. 
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VII. APPENDIX: CANONICAL TRANSFORMATION ON 
HUBBARD MODEL 

The goal of the perturbative approach is to bring 
the Hubbard model into block diagonal form in which 
each block has a fixed number of 'Active' doubly occu- 
pied sites. We say 'fictive' because the operators which 
make double occupancy a conserved quantity are not 
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FIG. 20 Carrier density computed from the spectral function 
shown in Fig. (|15|) . Clearly shown are two components: 1) 
one independent of temperature and scaling with the doping 
level and 2) the other temperature dependent describing the 
new composite or bound degrees of freedom. Fitting the car- 
rier density to Eq. ((7| enables an extraction of the gap A 
shown in the inset. The experimental values are also shown 
for LSCO: sol i d tria ngles (|Ando et a ll '2004bl :IOno a/.Ll2007l : 
iPadiUa et aH 120051 ) and sauaresljNishikawa aZ.LIl994l ') The 
excellent agreement obtained with the experimentally deter- 
mined values for the pseudogap indicates that the binding 
energy scale in the carrier density is the pseudogap energy. 



the physical electrons but rather a transformed (dressed) 
fer mion we cal l fjg d efined below. Following Eskes et 



al.( Eskes et al 



19941) ■ for any operator O, we define O 
such that O = 0(c) and O = 0(/), simply by replacing 
the Fermi operators Ciu with the transformed fermions 
ficr- Note that O and O are only equivalent in the 
U = oo limit. The procedure which makes the Hubbard 



model block diagonal i s now well known (|Eskes_ei 



11994 iMacDonald et all Il988l ). One constructs a simi 



larity transformation S which connects sectors that differ 
by at most one 'fictive' doubly occupied site such that 



H = e" He- 



rn 



becomes block diagonal, where H is expressed in terms 
of the transformed fermions. In the new basis, [H, V] — 
0, implying that double occupation of the transformed 
fermions is a good quantum number, and all of the eigen- 
states can be indexed as such. 

Our focus is on the relationship betwee n the physi- 
cal and 'fictive' fermions. To leading order (jEskes et all 



(35) 



are linear combinations of multiparticle states in the 
transformed basis as is expected in degenerate pertur- 
bation theory By inverting this relationship, we find that 
to leading order, the transformed operator is simply. 



fia 



(36) 



where 



X,, 



{tIju ^i(j^Cj(j Cj^CifjCifj -\- Cj^^CifjCja .(37) 



What we would like to know is what do the transformed 
fermions look like in the lowest energy sector. We accom- 
plish this by computing the projected operator 



(1 - nis)fia ~ (1 - n.iij)cia 4- — '^gij [(1 - nis-)X 



U 



(38) 



Simplifying, we find that 

(1 - ni^)f^a- ~ (1 - ni^)cig + ^Vac\gbi 



+ (1 - njff) [c^^gCia - CjacL) C^a -(39) 



-uj.ujja where j is 



Here V„ = -Vg ^ I and h, = J2ja^<yf^<^f: 
summed over the nearest neighbors of z. As is evident, 
the projected 'fictive' fermions involve the projected bare 
fermion, (1 — •nia)cia, which yields the 2x sum rule plus 
admixture with the doubly occupied sector mediated by 
the t/U corrections. These t/U terms, which are entirely 
local and hence cannot be treated at the mean-field level, 
generate the > 2x or the dynamical part of the spectral 
weight transfer. This physics (which ha s been shown to 



play a significant role even at half-filling ()Delannov et al. 



2OO5I )) is absent from projected i nodels such as the stan - 



dard implementation (|Andersonl 120061 : iLee et all 120061 ) 

of the t — J model in which double occupancy of bare 
electrons is prohibited. 

Consider now the low-energy Hamiltonian in the bare 
electron basis. The answer in the transformed basis 



is well-known (|Eskes et al\ . 119941 ) and involves the spin- 
exchange term as well as the three-site hopping term. 
Our interest is in what this model corresponds to in terms 
of the bare electron operators which do not preserve 
double occupancy. To accomplish this, we simply undo 
the similarity transformation after we have projected the 
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transformed theory onto the lowest energy sector. Hence, 



the quantity of interest is Hs 



Since in the transformed basis all such subspaces lie at 



least U above the m = sector, it is sufficient to focus 
on Pqc^ He~^ Pq. To express Poe^He~^Po in the bare 
electron operators, we subs titute Eg. into the first 



of Eqs. (14) of Eskes, et al. (|Eskes et aLl . ll994|) to obtain 



m 



e'^Poe^He^^Poe^ 



U ^ 



E 



{da [(1 - + + da^^<rV3^] + h-C.} (40) 



as the low-energy theory in terms of the original electron 
operators. Here, — Cia{l ~ riig-) and rji, 



VIII. APPENDIX: SPECTRAL FUNCTION AT FINITE 
DOPING 

Our analysis of the half-filled system points to a poten- 
tial organizing principle of the strong correlations. The 
dynamics leading to the Mott gap are completely inde- 
pendent of the spin-spin interaction contained in the term 
jfop. In fact, this term is strictly irrelevant relative to the 
terms of the form ip'^h. This suggests that we can drop 
the |&[^ term, without losing any relevant physics, as long 
as we are interested in the charge dynamics in the normal 
state. The resultant Lagrangian 



^Mott = E^"*^ ~ ni,-a)cl ,jCi^a- " * E ' 



f2 



h.c. , (41) 
I 



is quite simple in this case. In obtaining Eq. (|1T|) . we 
retained only the leading term in the t/U expansion for 
the fermionic matrix and of course dropped the |&p 
term from Eq. (PH). We refer to this Lagrangian as LMott 
as it contains s olely the charge dynamics. While the pre- 
dominant view! Anderson . 2006 : Lee et all . 12006 ) is that 
the spin-spin interaction dictates the physics in the Mott 
state, simple power counting in the exact low-energy the- 
ory reveals otherwise. Our analysis of the spectral func- 
tion of a doped Mott insulator is based entirely on Eq. 



We first assume that ip has no bare dynamics and is 
spatially homogeneous. The justification for this assump- 
tion is simple: 1) there are no gradient terms in the La- 
grangian for the charge 2e boson and 2) its primary role 
is to mix the sectors which differ in the number of dou- 
bly occupied sites. To proceed, we will organize the cal- 
culation of G{k,uj) by first integrating out the fermions 
(holding If fixed) 



Gik,Lu) = iy [Dip*][Dip]FT (^j [i?c*][i?c,]Q(t)c*(0)exp-/^"°"l^'^l'^* 



(42) 



To see what purely fermionic model underlies the 
neglect of the spin-spin term in Eq. we inte- 

grate over ipi in the partition function. The full de- 
tails of ho w to carry out such an integration are detailed 



elsewhere ljChov et all l2008bl ). The resultant Hamilto- 



nian is not the Hubbard model but rather. 



H' 



HHubb - ^l^i 



(43) 



a t — J — U model in which the spin-exchange interac- 
tion is not a free parameter but fixed to J —t^/U. 
That the t — J — U model with J — —t^/U is equivalent 
to a tractable IR model, namely Eq. (|4T|) without the 
spin-spin term, is an unexpected simplification. As Mott 
physics still pervades the t—J — U model in the vicinity 
of half- filling, our analysis should reveal the non-trivial 
charge dynamics of this model. The effective Lagrangian 
can be diagonalized and written as 



ka k kcr 



(44) 
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in terms of a set of Bogoliubov quasiparticles, 

7fct = + cos Okcl^ + sin OkC^ki (45) 
Jki = - sin Okcl^ + cos 9kC^ki (46) 

which define the propagating degrees of freedom in a hole- 
I 



doped Mott insulator, where cos^ = ■^{^ + Here, 

Qffe = 2(cosfc2,+cosfcy), Eq = -(2fi+^), Ek = -gttak- 

11 \k-- ^/El + Al Afe = s<fi*{l ~ fafe), and gt = Tfs 
when S— I— n^l — Q + 2(p*(p is a renormahzed factor 
which originates from the correlated hopping term (1 — 
nia-)cl^Cja{l — njg). Starting from Eq. (j44|) . we integrate 
over the fermions in Eq. (|42|) to obtain. 



G(k,io) = i / [i?(^*][i5<^]G(fc,^,^)exp-S.(^^o+£.-A,-|i„(i+cxp-^-'=)) (47) 



where 



G{k,uj,ip) = 



sm 



cos 



UJ + \k[ 



(48) 



is the exact Green function corresponding to the La- 
grangian, Eq. (|44l) . The two-pole structure of G{k,uj, ip) 
will figure prominently in the structure of the electron 
spectral function. To calculate G{k,Lu), we numerically 
evaluated the remaining ip integral in Eq. (H71) . Since 
Eq. (j47p is averaged over all values of (/?, we have circum- 
vented the problem inherent in mean-field or saddle-point 
analyzes. Physically, Eq. (|T7)) serves to mix (through the 
integration over ip) all subspaces with varying number of 
double occupancies into the low-energy theory. Hence, 
it should retain the full physics inherent in the bosonic 
degree of freedom. 
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